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psj . We investigate analytically and numerically the role of quantum fluctuations in 

I reconstruction of optical objects from diffraction-limited images. Taking as exam- 
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pie of an input object two closely spaced Gaussian peaks we demonstrate that one 
can improve the resolution in the reconstructed object over the classical Rayleigh 
limit. We show that the ultimate quantum limit of resolution in such reconstruction 
procedure is determined not by diffraction but by the signal-to-noise ratio in the 
' input object. We formulate a quantitative measure of super-resolution in terms of 
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T^lj- , the optical point-spread function of the system. 

o 



PACS numbers: 42.50.Dv, 42.30.Wb, 42.50.Lc 



§ ■ I. INTRODUCTION 

Quantum imaging is a new branch of quantum optics that investigates the ultimate 



X 

' performances of optical imaging allowed by the laws of quantum mechanics 



IS- One 

of the pertinent questions of quantum imaging is about the quantum limits of the optical 
resolution. The classical resolution limit was established by Abbe and Rayleigh at the end 
of the nineteen century and is very well known nowadays. The Rayleigh resolution criterion 
states that the resolution in an optical system is limited by diffraction on its optical elements 
due to the wave nature of the light. According to this criterion, two closely spaced points at 
the input of the optical system cannot be resolved if the distance between them is smaller 
than A/NA where A is the wavelength of the light and NA - the numerical aperture of the 
optical system. 

This classical resolution limit is based on the presumed resolving capabilities of the human 
eye and is not a fundamental limit like, for example, the Heisenberg uncertainty relation. 
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Nowadays using modern CCD cameras for detection of optical images with subsequent elec- 
tronic processing one can often improve the resolution beyond the classical diffraction limit. 
For example, while the Rayleigh criterion puts the resolution limit of about 0.2 fim for an 
optical microscope, by processing microscopic images one can achieve a precision of about 
Inm jSl. Even more spectacular results has been obtained recently for detection of a small 
displacement of a laser beam Using "spatially squeezed" laser beam the authors of |^ 
have succeeded to measure a transverse displacement of 1.6 Angstrom of a laser beam with 
wavelength A = 1064 nm. 

These two experiments are particular examples of the so-called super-resolution tech- 
niques that aim to improve the optical resolution beyond the classical diffraction limit. 
Super-resolution is possible when one has some a priori information about the object. For 
example, in practice one usually deals with objects of finite size. In this case one can ob- 
tain super-resolution because the spatial Fourier spectrum of such objects, produced in the 
focal plane of a lens where the system pupil is located, is an analytical function. Thus, 
measuring a part of the Fourier spectrum, transmitted by the system pupil, one can in 
principle reconstruct the rest of the spectrum via the analytical continuation and, therefore, 
obtain an infinite resolution. However, the procedure of such an analytical continuation is 
extremely sensitive to different kind of noises present in the optical system. Recently it was 
shown p that the ultimate limit of resolution in diffraction-limited imaging is determined 
not by diffraction but by the quantum fluctuations of light within the object area and the 
vacuum fluctuations outside it. These quantum fluctuations set up the standard quantum 
limit of resolution which can be by the orders of magnitude smaller than the diffraction 
limit. Moreover, one can further improve the resolution beyond the standard quantum limit 
using multimode squeezed light. An optical scheme of such a super-resolution microscopy 
with squeezed light was proposed in Ref. ^j. 

In this paper we numerically simulate the role of quantum fluctuations in reconstruction 
of optical objects from the diffraction-limited images. We take as example of an input object 
two Gaussian peaks placed so close to each other that in the output image they cannot be 
resolved according to the Rayleigh criterion. We assume that one can measure the spatial 
Fourier components of the input object in the focal plane of the imaging lens within the pupil 
area. This kind of measurement corresponds to the super-resolving Fourier-microscopy [7| 
since instead of observing the image one observes its spatial Fourier spectrum. We formu- 
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late the quantum theory of this microscopy in terms of the prolate spheroidal functions 
very similar to the theory developed in Ref. Q]. Using this quantum theory, we simulate 
numerically the reconstructed objects and investigate the role of quantum fluctuations on 
their resolution. 

Our numerical simulations allow us to confirm that when quantum fluctuations are not 
taken into account, one can easily improve the resolution in the reconstructed objects over 
the diffraction limit by one order of magnitude. When the object is illuminated by a light 
wave in coherent state, the quantum fluctuations of light inside the object and vacuum 
fluctuations outside it set up the standard quantum limit of resolution which depends on 
the signal-to-noise ratio in the input object. Finally, we demonstrate that one can go beyond 
the standard quantum limit of resolution using the multimode squeezed light for illumination 
of the object. 

The paper is organized as follows. In Section |n] we formulate the quantum theory of 
the super-resolving Fourier-microscopy in terms of prolate spheroidal wave functions. In 
Section IIIII we give numerical examples of super-resolution of the input object consisting 
of two closely spaced Gaussian peaks and illustrate qualitatively the role of quantum fluc- 
tuations on the degree of super-resolution. In Section |TV| we provide quantitative measure 
of super-resolution in terms of the optical point-spread function. In particular we plot 
super- resolution factor as a function of the mean photon number in the input object for 
the coherent light and multimode squeezed light. We discuss our results and the future 
perspectives in Section IVl 

II. QUANTUM THEORY OF THE SUPER-RESOLVING 
FOURIER-MICROSCOPY 

In this section we shall present the original quantum theory of the super-resolving Fourier- 
microscopy in terms of prolate spheroidal functions. The optical scheme of diffraction-limited 
coherent optical imaging is shown in Fig. Q For simplicity we consider one-dimensional case. 
The object of finite size X is placed in the object plane. The first lens Li performs the spatial 
Fourier transform of the object into the pupil plane with a pupil of finite size d. Diffraction 
on this pupil is a physical origin of the finite resolution in our scheme (we neglect diffraction 
on the imaging lenses). The second lens L2 performs the inverse Fourier transform and 
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creates a diffraction-limited image in the image plane. 

As mentioned above, to achieve super-resolution one needs some a priori information 
about the object. In our case we know a priori that the object is confined within the area of 
size X and is identically zero outside. The spatial Fourier transform of such an object is an 
entire analytical function. Therefore, knowing the part of the Fourier spectrum within the 
area d of the pupil allows for an analytical continuation of the total spectrum and, therefore, 
for unlimited resolution. However, this analytical continuation is extremely sensitive to the 
noise in the diffracted image and, as shown in Ref. j^, is limited by the quantum fluctuations 
of light. 

Quantum theory of the diffraction-limited optical imaging in Fig. ^ was developed in 
Ref. 0] in terms of prolate spheroidal wave functions. These are the eigenfunctions of the 
imaging operator describing the transformation of the optical fleld from the object plane 
into the image plane. One can decompose the input object and the output image over 
these eigenfunctions and obtain the relation between the decomposition coeflicients. Then 
detecting the output image with, for example, a sensitive CCD camera, one can evaluate 
the decomposition coefficients of the image. Using the relation between the decomposition 
coefficients of the image and the object one can reconstruct the latter with resolution better 
than the classical diffraction limit. 

Our numerical simulations in Ref. have shown, however, that for evaluation of the 
decomposition coefficients one has to detect the output image over unrealistically large area 
in the image plane due to the oscillating behavior of the prolate functions. That is why in 
Ref. we have proposed a modified version of the scheme where the CCD camera is placed 
in the pupil plane instead of the image plane, and one detects the spatial Fourier spectrum. 
We have called this modified scheme super-resolving Fourier-microscopy. The advantage 
of this set-up is that now the spatial Fourier spectrum is measured over the finite region 
within the pupil. To understand the role of the quantum fiuctuations on the resolution of 
the Fourier-microscopy we need to formulate the quantum theory of this modified scheme. 

Let us introduce the dimensionless spatial coordinates in the object plane as s = 2x/X, 
and in the pupil plane as = 2y/d (see Fig. 1). The dimensionless photon annihilation 
operators in the object plane will be denoted as a(s) and in the pupil plane as /(O- These 
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operators obey the standard commutation relations, 

[a(.), at(.')] = 5{s - s'), [/(O, PiC)] = ^(^ - (2-1) 

and are normalized so that (d^(s)d(s)) gives the mean photon number per unit dimensionless 
length in the object plane and " pupil plane. The spatial Fourier transform 

{Td){C,) performed by the lens Li, in terms of these dimensionless variables reads as follows, 



m) = (Tarn = / a{s)e—^ds, (2.2) 



7r 

where c = — -ry is the space-bandwidth product of the imaging system. It is important to 
2 A/ 

note that the limits of integration in this equation are over the whole object plane since it 



concerns operators and not the classical c-numbers. 

As in Ref. jJl we shall formulate our quantum theory of super-resolving Fourier- 
microscopy in terms of prolate spheroidal functions ipkis) These are the eigenfunction 
of the imaging operator of the scheme, orthonormal on the interval —oo<s<oo. To ob- 
tain the canonical transformation of the photon annihilation and creation operators from 
the object into the pupil plane we shall split the coordinates s and ^ into two regions, the 
"core", |s| < 1 and |^| < 1, corresponding to the area of localization of the classical object 
and the transmission area of the pupil, and the "wings", |s| > 1 and |^| > 1, outside these 
areas. The orthonormal bases in these areas of the object plane are given by 
1 



ipkis) = { Xk{s) = { 1 , _ ■ ■ ^ (2.3) 



Ms) \s\ < 1, 1^1 < 1' 



1^1 >!' [VT^ 



ipkis) \s\ > 1, 



k 



where are the eigenvalues of the corresponding prolate spheroidal functions ipkis), de- 
pending on the space-bandwidth product c. We should note that the functions (fk{s) are 
complete in the Hilbert space L^(— 1, 1). Similar relations take place in the pupil plane. 

In terms of two sets and {Xk{s)} we can write the annihilation operators in the 

object plane as 



and in the pupil plane as 



a{s) = Y.^kVk{s) + ^bkXk{s), (2.4) 

k=0 k=0 



ho = E hVkiO + E dkXkiO- (2.5) 

k=0 k=0 
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Here and fk are the annihilation operators of the prolate modes (pk in the core region of 
the object and the pupil planes, while bk and cjk are the annihilation operators of the prolate 
modes Xk in the wings regions. The operators dfc and bk are expressed through the field 
operator d(s) by 



ak 



a{s)(pk{s)ds, 



a{s)xk{s)ds. 



(2.6) 



Similar relations hold for fk and gk and /(O- 

Four sets of operators dk, bk, fk, and gk obey the standard commutation relations of 
the photon annihilation and creation operators of the discrete modes. For example, for the 
operators dk in the core region of the object plane we have. 



dk,dj^,] — 6kk', [dk,dk'] — 0. 



(2.7) 



The same relations are fulfilled for three other groups of operators bk, fk, and cjk. It is clear 
that creation and annihilation operators of different groups are independent and commute 
with each other. 

In our analysis we shall use the following properties of prolate spheroidal functions [lol |. 



Ms)e-"''^ds 



(2.8) 



Ms)e-'''^ds = {-t)^J—cpkiO, (2.9) 

V c 

Using (j2.3p . the field transform (j2.2p between the object and the pupil plane and these 
properties, we find the following propagation relations for the core and wings of the light 
wave, 



{TXk) (0 = 



(2.10) 
(2.11) 



1 - AfcV5fc(s) - \/AfcXfc(s) 

Substituting Eqs. ()2.4|) and ()2.5p into the field transform ()2.2|) and using Eqs. ()2.10|) and 
(12.1111 , we arrive at the following relations between the photon annihilation operators of the 
prolate modes in the object and the pupil planes. 



fk = {-i) ( V AfeSfc + V 1 - Xkbk) 



9k = {-i) (yl - AfcOfc - J Xkbk)- 



(2.12) 
(2.13) 
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This relations are similar to the transformation performed by a beam-splitter with the 
amplitude transmission coefficients {—i)''\/Xk and the reflection coefficients — Xk, 

and preserves the commutation relation of the annihilation and creation operators in the 
pupil plane. 

Let us assume that we can detect the spatial Fourier amplitudes f{^) in the pupil plane 
within the transmission area of the pupil using a sensitive CCD camera. This transmitted 
part of the spatial Fourier spectrum is given by the first sum in Eq. ()2.5|) . the term given 
by the second sum is absorbed by the opaque area of the pupil. It should be emphasized 
that, since we need the complex field amplitudes and not the intensities, one should use the 
homodyne detection scheme with a local oscillator. Using Eqs. ()2.5j) . ()2.(j|l and ()2.12j) we can 
calculate the operator-valued coefficients a^*^^ of the reconstructed object as 

where the superscript (r) stands for "reconstructed". As follows from Eq. ()2.14j) . the recon- 
struction of the input object is not exact because of the second term in Eq. ()2.14|) . This 
term contains the annihilation operators responsible for the vacuum fluctuations of the 
electromagnetic fleld in the area outside the object. It is important to notice that these 
vacuum fluctuations prevent from reconstruction of the higher and higher coefficients dfc in 
the object because of the multiplicative factor ^^(1 — Xk)/Xk- Indeed, the eigenvalues 
become rapidly very small after the index k has attained some critical value. This leads to 
rapid "amplification" of the vacuum fluctuations in the reconstructed object that limits the 
number of the reconstructed coefficients dk- 

III. QUANTUM FLUCTUATIONS AND RECONSTRUCTION OF THE 
SPATIAL FOURIER SPECTRUM OF THE OBJECT 

A. Reconstruction of classical noise-free objects 

In this section we shall illustrate numerically the role of quantum fluctuations on the 
reconstruction of simple objects with super-resolution beyond the classical diffraction limit. 
However, before taking into account quantum fluctuations of light in the input object, 
we would like to demonstrate the potential of the super-resolution technique with prolate 
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spheroidal functions for reconstruction of noise-free classical objects, i. e. when the quan- 
tum fluctuations are neglected. This case corresponds to the classical limit of the quantum 
theory developed in previous section and can be simply obtained by taking mean values of 
the operators. 

In what follows we shall denote the classical complex amplitudes corresponding to the 
quantum-mechanical operators by the same letters without carets, for example a{s) = {a{s)), 
flfc = {O'k), et cetera. Since the classical complex amplitude of the object is zero outside the 
area |s| < 1, we have (bk) = 0. Using Eq. ()2.4|) we can write this classical amplitude as 

oo 

= E «fc<^fc(s). (3.1) 

The classical complex amplitude /(^) of the field in the pupil plane is obtained from Eq. ()2.2|1 . 

^(^) = \/? ais)e-^^'^ds, (3.2) 
V 2tt j-1 

with the integration limits over the object area, |s| < 1. Taking into account the property 
of the prolate functions given by Eq. ()2.8|) . we can write the spatial Fourier spectrum /(^) 
as the following decomposition, 

oo 

/(O = E(-0Wfc(5), -oo<^<oo. (3.3) 

k=0 

This spatial Fourier spectrum of the object spreads outside the transmission area of the 
pupil 1^1 < 1. The spatial Fourier components in the opaque area are absorbed and cannot 
be detected by the CCD camera placed in the pupil plane. Super-resolution attempts to 
reconstruct these absorbed Fourier components. From Eq. ()2.14|1 we obtain the classical 
reconstructed coefficients, 

a^^ = ak. (3.4) 

Because we have neglected the quantum fluctuations, the reconstructed coefficients are iden- 
tical to those of the input object. The classical amplitude of the reconstructed object a'^''"^(s) 
can be written as the following decomposition over the prolate functions, 

L-l 

a('")(s) = ak(pk{s). (3.5) 

fc=0 

Since in practice one can never have infinitely many coefficients a^, we have restricted the 
summation in this equation to L first prolate functions. When L ^ oo, the reconstructed 
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object approaches the exact one, a^^^(s) — > a{s). In practice the super-resolution over the 
Rayleigh hmit is determined by the number L of terms used in the decomposition ()3.5|1 . 

Alternatively to reconstruction of the object itself one can try to reconstruct its spatial 
Fourier spectrum as 

f^^\0 = j:i-')'^kMs), -oo < e < oo. (3.6) 

fc=0 

Similarly to the reconstruction of the object, when L oo, the reconstructed spectrum 
approaches the exact one, /'■^-'(O ~^ /(O- 

For numerical simulations we have taken a simple object of two narrow Gaussian peaks, 

aW = A[exp(-(i^)+exp(-(i±^)]. |.| < 1. (37) 

of width a separated by distance 2so. We choose 2sq = 1 and a = 0.1, so that two peaks 
are well separated in the input object. The normalization constant A is chosen so that the 
integral of the object intensity over the area of the object is equal to the total mean number 
of photons {N) in the object, 

J^^a\s)ds = (N). (3.8) 

The Rayleigh resolution distance R = ttX/ (2c) in dimensionless coordinates s is equal to 
tt/c, where c is the space-bandwidth product. In our simulations we work with c = 1. In 
this situation for 2so < tt we are beyond the Rayleigh limit. 



In Fig. 2a we have plotted the normalized input object a{s) / y (N), in Fig. 2b - its spatial 
Fourier spectrum f{^) in the pupil plane, and in Fig. 2c - the output image e{s) in the image 
plane. Comparing the input object with its image one can clearly see that it is impossible 
to resolve two Gaussian peaks in the image plane according to the Rayleigh criterion. In 
Fig. 2b we have shown by the grey color the opaque area of the pupil. The part of the spatial 
Fourier spectrum of the object in this area is absorbed and therefore cannot be detected by 
the CCD camera placed in the pupil plane. Below we shall illustrate the reconstruction of 
these absorbed Fourier components by the technique of the prolate functions. 

For numerical simulations we had to evaluate two sets of prolate spheroidal functions, 
<fk{s), defined on the interval |s| < 1, and ipkis) defined for all s, — oo < s < oo. The first set 
is necessary for decomposition of the input object a(s), while the second one is needed for 
reconstruction of the spatial Fourier spectrum /'•'"■' (0. For numerical calculations of ^k{s) 
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we have used the algorithm from Ref. ^|. In this algorithm the prolate functions fk{s) are 
evaluated as the series with the Legendre polynomials Pk{s), 



OO ( 

Ms) = T.lt\k+-Pk{s), \s\<l. (3.9) 
k=o * ^ 

The coefficients 7^"'' are found as the eigenvectors of the symmetric matrix A with the 



following nonzero elements, 



, (k + 2)(k+l] 



A.k,k+2 — Ak+2.k — C , , (3-11) 

(2A; + 3)J(2A; + l)(2A; + 5) 



for all A; = 0,1,2. . . 

We have written a numerical program in MATHEMATICA which implements this algo- 
rithm. The advantage of this method is that it does not require the direct solution of the 
eigenproblem for ipk{s) and which is unstable due to the rapid decrease of the eigenval- 
ues. This algorithm allows us to calculate for c = 1 at least 17 first prolate functions in 
spite of the fact that the eigenvalues of the higher order functions become extremely small 
(for example, A17 = 4.183 x 10~^°). In Fig. 3 we show the first 17 prolate functions (pk{s) 
evaluated by our numerical program. 

For numerical calculation of the second set of the prolate spheroidal functions, il>k(s), we 
have used the following property of the Legendre polynomials ( see Eq. (10.1.14) in Ref. (l2j|). 

£^ Pn{s)e-'''^ds = 2z"j„(0, (3.12) 

where j„(s) is the spherical Bessel function of the first order Using this equation we 
can easily obtain the following representation of ipn{C}-> 



MO = ^"f:(-0'7i"Y^+ -00 < e < 00. (3.13) 

We illustrate the result of reconstruction of the spatial Fourier spectrum of the input object 
in Fig. 4. In this figure we show the exact spatial Fourier spectrum of the input object, 
drown by a solid line, as a function of dimensionless coordinate ^ in the pupil plane. Only 
part of this spectrum within the transmission area of the pupil, |.^| < 1, is transmitted to 
the image plane. The spatial Fourier harmonics in the opaque area of the pupils, shown by 
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grey color, are absorbed. This is a reason of very large diffraction spread in the image plane 
shown in Fig. 2c. The dashed lines in Fig. 4 correspond to the spatial Fourier spectrum 
of the reconstructed object with L = 5,7 and 11 prolate functions. One can see that the 
reconstructed spectrum approaches the exact one for ever higher spatial frequencies |^| as 
the number of prolate functions increases. 

With 7 prolate functions two spectra are very close to each other for spatial frequencies 
\^\ < 8. This corresponds to a super-resolution factor of 8 over the Rayleigh limit. 

B. Reconstruction of objects with quantum fluctuations 

For numerical simulations of quantum fluctuations we have chosen a c-number represen- 
tation of the quantum mechanical operators and bk in Eq. ()2.14|) corresponding to the 
antinormal ordering of the creation and annihilation operators. In this representation the op- 
erators dfc and bk become the c-number Gaussian stochastic variables ak and (3k respectively, 
which we shall write as 

ak = ak + Sak, (3k = S(3k. (3-14) 

Here a^, = {hk) is the mean value of the field coefficients in the object area, and 5ak and 5(3k 
are the stochastic Gaussian fluctuations. Note that the mean values {bk) are zero because 
the classical field component outside the object vanish. We have chosen the antinormally- 
ordered representation because it remains valid even in the case of the multimode squeezed 
state of light field at the input of the scheme. 

We introduce the quadrature components of the fiuctuations 5ak and 513k as follows, 

5ak = 6X^ + i5Y^, 6 (3k = 5Xf + i6Y^^. (3.15) 

When the input light is in the coherent state in the object area and in the vacuum state 
outside, the correlation functions of the quadrature fiuctuations are equal to 

(SX^SX^,,) = {6Y,^5Y^:) = ^Skk', (3.16) 

with fi = a, p. 

If instead of coherent light we use the multimode squeezed light for illumination of the 
object and multimode squeezed vacuum in the area outside with subsequent homodyne 
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detection at the pupil plane, these correlation functions become 

{5Xi:5Xj:,) = {5Y,'5Y,^,) = (3.17) 

where r is the squeezing parameter. In these formulas we have assumed that the input light 
is amplitude-squeezed and for simplicity have chosen the same squeezing parameter for all 
the essential modes that are used in the decomposition of the reconstructed object. 

The relative value of quantum fluctuations depends on the signal-to-noise ratio in the 
input object which for the light in a coherent state is determined by the total mean number 
of photons passed through the object area during the observation time. For example, for a 
laser beam with A = 1064 nm and optical power of 1 mW, and observation time of 1 ms we 
obtain the mean photon number of (N) = 5.3 ■ 10^^. 

In Fig. 5a we have shown the results of reconstruction of the spatial Fourier spectrum 
of the object from Fig. when quantum fluctuations of a coherent state are taken into 
account. The solid line gives the exact spatial Fourier spectrum of the object. As in Fig. 4 
the grey area shows the absorbing part of the pupil. We use 7 prolate functions and the 
mean photon number in the input object is taken (N) = 10^^. The five thin lines correspond 
to the five random Gaussian realizations of the quantum fluctuations in the coherent state 
of dk and the vacuum fluctuations of bk- The dashed line corresponds to the reconstructed 
spectrum with 7 prolate functions without noise. One can observe that the role of quantum 
fluctuations becomes more and more important as one goes to the higher and higher spatial 
frequencies where the random realizations of the Fourier spectra deviate more and more 
from the mean value given by the dashed line. 

In Fig. 5b we have increased the total mean value of photons to (A^) = 10^'^. This 
corresponds to an increased signal-to-noise ratio in the input object and should allow for 
better super-resolution. This is illustrated in Fig. 5b by reduced deviation of the random 
realizations from the mean value of the spectrum a compared to Fig. 5a. 

The same result can be achieved by using multimode squeezed light instead of increasing 
the power of the source illuminating the object. This is illustrated in Fig. 5c where we have 
used (A^) = 10^^ as in the Fig. 5a, but have considered the light in a multimode squeezed 
state with the squeezing parameter e'" = 10 instead of the coherent state. As the result the 
fluctuations in the higher spatial frequencies are reduced that gives better super-resolution. 

In next section we shall give a quantitative characteristic of super-resolution as a function 
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of the signal-to-noise ratio. 

IV. POINT-SPREAD FUNCTION FOR SUPER-RESOLVING 
RECONSTRUCTION OF OBJECTS 

In modern classical optics the resolution of an optical system is characterized not by the 
two-point Rayleigh resolution criterion, but in terms of its spatial transmission bandwidth. 
A typical optical system has a finite band of spatial frequencies that are transmitted through 
the system up to some cut-off frequency determined by the size of the system pupil. The 
optical system is then said to be bandlimited or diffraction-limited since diffraction effects 
on its pupil are responsible for finite resolution. 

A coherent diffraction-limited imaging system in classical optics can be described by a 
linear equation relating the complex amplitude a{s) of an input object with the complex 
amplitude e{s) of the image 



The impulse response function h{s, s') that appears in this integral equation represents the 
image at point s in the image plane from a point-source at point s' in the object plane. For 
translationally invariant or isoplanatic systems the impulse response depends only on the 
difference s — s' and the integral in ()4.H1 becomes convolution, 



In optics, the impulse response h{s — s') is usually called the point-spread function (PSF) 
of the system, and its Fourier transform the transfer function (TF). For bandlimited optical 
systems the transfer function is identically zero outside the transmission band of the system. 
Super-resolution is defined as technique of restoring the spatial frequencies of the object 
outside the transmission band It is important to underline that in case when the object 
and the image fields are related by the convolution ()4.2|1 . super- resolution is impossible. To 
achieve super-resolution one needs some a priori information about the input object. In our 
case the a priori information is the assumption that the object has finite size. With this 
assumption Eq. ()4.2p in dimensionless coordinates becomes jsj. 




(4.1) 




(4.2) 




(4.3) 
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with the imaging PSF h{s — s') given by 

h(s - s') = — ^-^ (4.4) 

For the reconstruction process we can write similar relation between the reconstructed field 
operator a^^\s) and the object field operator d(s). Using an operator- valued equivalent of 
Eq. ()3.5j) together with Eq. ()2.14j) we arrive at the following result, 



aW(s) = /' h^^\s,s')d{s')ds' +J2 J^-^hMs). (4.5) 
Here the reconstruction point-spread function h^'^\s,s') is given by 

h^-\s,s') = J2Ms)Ms')- (4.6) 

fc=0 

As seen from this equation, the form of the reconstruction PSF and, in particularly, its width 
depends on the number of terms L in the sum. When this number grows infinitely, L ^ oo, 
the reconstruction PSF tends to 5-function, 

oo 

hm /iM(s, s') = J2 Ms)Ms') = S{s - s'), (4.7) 

and we have unlimited super-resolution. However, this ideal situation is never realized 
practically due to the second term in Eq. ()4.5j) which grows infinitely when L ^ oo. Thus, 
Eq. ()4.5|) is a good illustration of the statement that the ultimate limit of super-resolution 
in the reconstructed object is given not by diffraction but by the quantum fluctuations of 
light represented by the second term. 

The number L of terms in the sum ()4.fij) which determines the width of the reconstruction 
PSF, depends on the signal-to-noise ratio in the input object. To obtain the maximum L we 
shall compare the signal-to-noise ratio in the input object to that in the reconstructed object. 
As follows from Eq. ()4.5|) with increasing L the signal-to-noise ratio in the reconstructed 
object deteriorates. We shall assume that reconstruction of the object is possible until the 
limit when the signal-to-noise ratio in the reconstructed object becomes unity. 

Let us deflne the singal-to-noise ratio in the input object as 3| 

R^^, (4.8) 

{(ANr-) 

where 

(AT) = y\a^(s)a(s))ds, (4.9) 
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is the total mean number of photons in the input object, and ((AA^)^) - its variance. Similarly 
we define the signal-to-noise ratio R^"^^ in the reconstructed object as 

i?M = ^ . ' — , (4.10) 



where the mean number of photons in the reconstructed object is given by 

{N^'^) = l^^{a^'^^{s)a^'\s))ds. (4.11) 

The deterioration of the signal-to-noise ratio in the reconstructed object can be described 
by the noise figure F, 

F=^, (4.12) 

that is commonly used in the literature about amplifiers. Because the signal-to-noise ratio 
R^"^^ in the reconstructed object is always smaller than that in the input object, the noise 
figure is always larger than unity. If we assume that the minimum value of R^^^ that allows 
for reconstruction of the object is unity, this gives us the maximum noise figure -Fmax = R 
corresponding to the maximum super-resolution. 

Let us consider an input object in a coherent state, so that (a(s)) = a(s), (d"'"(s)a(s)) = 
|a(s)p, (a"''(s)a^(s')a(s')a(s)) = |a(s)p|a(s')p. It is easy to show that in this case the input 
signal-to-noise ratio R is equal to the mean total photon number in the input object, 

R= {N) = j^^ \a{s)\^ds. (4.13) 

On the other hand, for the signal-to-noise ratio i?'-^-' in the reconstructed object in this case 
we obtain the following result 



fc=0 /c=0 

where are the coefficients of decomposition of a^^\s) over the prolate functions 99^(5) in 
Eq. ((S3D. 

As follows from Eq. ()4.14j) . the signal-to- noise ratio and, therefore, the noise figure 
F depend on the shape of the input object. For numerical evaluation of the super-resolution 
factor as a function of the total mean number of photons in the input object we have taken 
a narrow rectangular object placed at the origin s = 0, 



UN) 

a.{s) = { (4.15) 
|s|>e/2. 
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Taking the width e of this object ever smaller we arrive at a point-like source, while keeping 
the total number of photons constant and equal to (N). Such a point-like object gives us 
the reconstruction PSF h^^\0, s) at the output. 

The degree of super-resolution in the reconstructed object can be characterized by the 
ratio of the width of the diffraction-limited imaging PSF to the width of the reconstruction 
PSF. In Fig. 6 we have shown the imaging PSF h{s) and the reconstruction PSF h^'^\0,s) 
ior L — 7 normalized to unity at their maxima. To define the super-resolution factor we 
shall introduce the half-widths W and Wl of these two PSF measured at their half-maxima. 
Then we define the super- resolution factor S as the ratio of W to Wl, 




(4.16) 



For the example given in Fig. 6 these half-widths are equal toW — 1.895, Wl — 0.252, and 
5 = 7.5. 

In Fig. 7 we have plotted the super-resolution factor 5" as a function of the total mean 
number of photons (N) in the input object for the case of coherent light and multimode 
squeezed light. As seen from this figure, for the same mean number of photons multimode 
squeezed light provides higher super- resolution that the coherent light. 



V. DISCUSSION AND PROSPECTIVES 



In this paper we have investigated analytically and numerically the quantum limits of 
super-resolution in reconstruction of optical objects from the diffraction-limited images. We 
assume that such a reconstruction is performed electronically from the data collected in a 
homodyne detection of images by means of a CCD camera placed either in the image plane 
or in the Fourier plane. We call the first situation super-resolving microscopy and the second 
- super-resolving Fourier-microscopy. In Section II we have presented the quantum theory 
of the latter. 

From our numerical simulations in Sections III and IV we conclude that a priori informa- 
tion about the finite size of the input object allows one to obtain significant super-resolution 
in the reconstructed object over the classical Rayleigh limit. For classical noise-free ob- 
jects we have demonstrated a possibility to achieve super-resolution of factor 8 over the 
Rayleigh limit. When quantum fluctuations in the object are taken into account, the degree 
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of super-resolution depends on the signal-to-noise ratio in the input object. When the ob- 
ject is illuminated by the light in a coherent state, this signal-to-noise ratio is given by the 
mean total number of photons during the observation time. We have quantitatively defined 
the super-resolution factor as a ratio of widths of the imaging and reconstruction point- 
spread functions. We have numerically evaluated this super-resolution factor as a function 
of the mean total number of photons for the coherent light and multimode squeezed light. 
For the same mean number of photons multimode squeezed light allows to achieve higher 
super- resolution than the coherent light. 

As future prospectives for further development of the theory presented in this paper we 
would like to mention the problem of read-out of binary information from optical discs, like 
CD and DVD, or optical memory. The storage density for the optical discs currently is 
limited by the spot size of the diffraction-limited focussed light beam. One possibility of 
increasing the storage density would be an attempt to put several bits of information inside 
the diffraction-limited light spot with subsequent super-resolution in the read-out of this 
information. In this case the amount of a priori information is clearly superior to the case 
considered in present article since for finite number of bits inside the laser spot one has a 
finite possible combination of light patterns. Therefore, one would expect higher potential 
of super-resolution in read-out of binary information as compared to the situation that we 
have investigated here. In the context of the problem of read-out of optical discs one would 
have to generalize our theory to optical systems with high numerical apertures which is 
usually the case for the systems of the optical data storage. 

This work was supported by the Project QUANTIM (IST-200-26019) of the European 
Union. 
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FIG. 1: Optical scheme of one-dimensional coherent diffraction-limited imaging. 
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FIG. 2: Double-peak object a(s) used in numerical simulations (a), its spatial Fourier spectrum 
/(^) in the pupil plane (b), and the image e(s) created in the image plane (c). Grey area shows 
the part of the spatial spectrum absorbed by the opaque area of the pupil. 
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FIG. 3: Examples of the prolate spheroidal functions <^fc(s) and the corresponding eigenvalues 
calculated using our numerical program. 
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FIG. 4: Exact spatial Fourier spectrum of the object from Fig. 2a (solid line), and the spectra 
reconstructed with L = 5, 7 and 11 prolate functions (dashed lines). Here and in Fig. 5 grey area 
indicates the absorbed part of the spatial spectrum. 
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FIG. 5: Exact spatial Fourier spectrum of the object from Fig. 2a (solid line), the reconstructed 

spectrum with L = 7 prolate functions (dashed lines), and five random Gaussian realizations of 
the reconstructed spectrum with L = 7 prolate functions (thin lines); (a) coherent light with mean 
total photon number (N) = 10^^, (b) coherent light with (N) = 10^^, and (c) squeezed light with 
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FIG. 6: DifFraction-limited imaging point-spread function ^(s) and the reconstruction point-spread 
function /i(^)(0,s) using 7 prolate functions. 




FIG. 7: Super-resolution factor 5 as a function of the total mean number of photons (TV) for 
coherent light (solid line, and multimode squeezed line with exp(r) = 10 (dashed line). 



